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Abstract. Near threshold neutral pion electroproduction on the deuteron is studied in the framework of 
heavy baryon chiral perturbation theory. We include the next -to-leading order corrections to the three- 
body contributions. We find an improved description of the total and differential cross section data mea- 
sured at MAMI. We also obtain more precise values for the threshold S-wave multipoles. We discuss in 
detail the theoretical uncertainties of the calculation. 
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1 Introduction 

Pion electroproduction off deuterium, j*(k) + d — > 7r°(g) + 
d, allows to extract the elusive elementary neutron ampli- 
tude. First data at small photon virtuality k 2 = —0.1 GeV 2 
and small values of the pion excess energy AW (AW = 
0.5, . . . , 3.5 MeV) have been taken at the Mainz Mikrotron 
MAMI-II [1]. In the threshold region, this reaction can 
be analyzed in the framework of chiral perturbation the- 
ory [2,3]. In Ref. [3], a (partial) next-to- leading order cal- 
culation was presented, which led to a satisfactory descrip- 
tion of the differential and total cross section data. That 
analysis can be improved in two respects. First, the so- 
called three-body corrections (meson exchange currents) 
were only considered at third order in that paper. Second, 
a similar remark holds for the single scattering proton and 
neutron P-waves. In this paper, we are concerned with the 
first issue, namely the next-to-leading (fourth) order cor- 
rections to the three-body corrections. As will be shown, 
this introduces (in principle) no new parameters and with 
that extension, the calculation is done with the same the- 
oretical precision as the fairly successful analysis of co- 
herent neutral pion photoproduction off the deuteron [4] . 
It also leads to an improved description of the data as 
shown in this paper. The second extension can not be 
done so straightforwardly. This is because so far only the 
fourth order corrections to the photoproduction P-wave 
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multipoles are available in the literature [5]. For the elec- 
troproduction case, similar calculations are underway [6], 
triggered in particular by the rather unexpected results for 
neutral pion electroproduction off the proton at the pho- 
ton virtuality k 2 = —0.05 GeV 2 [7]. These data can not be 
described within chiral perturbation theory at one-loop 
accuracy. This is in stunning contrast to the fairly good 
description of the data at higher photon virtuality from 
NIKHEF [8] and MAMI [9] that was obtained in [10]. Note 
that remeasurements at the low photon virtuality seem to 
be more in line with expectations from chiral perturbation 
theory or sophisticated models [11]. As we will argue, the 
effect of such single scattering P-wave contributions can be 
effectively included in a refit of some four-nucleon-photon 
operators, which leads to a much improved description of 
the data. Given the precision of the data, a more elaborate 
treatment appears inappropriate. 

The manuscript is organized as follows. In Section [2 
we briefly discuss the pertinent formalism. We heavily bor- 
row from Ref. [3] and only discuss the new fourth order 
three-body contributions in some detail. Our results are 
presented and discussed in Section |31 and we end with a 
short summary and outlook in Section 0] Some technical- 
ities are relegated to the appendices. 



2 Effective field theory description 

2.1 General remarks 

In Ref. [3] , we developed the multipole formalism for near 
threshold pion electroproduction off deuterium (see also 
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[12,13]) and calculated the pertinent transition matrix el- 
ements in chiral perturbation theory. This is based on an 
effective chiral Lagrangian of pions and nucleons chirally 
coupled to external sources like the photon, 

£cff = -Ctttt + £irN + £nn , (1) 

where the first term subsumes the interactions between 
the Goldstone bosons, the second the pion-nucleon inter- 
actions and the third term the short range part of the 
two-nucleon interaction. All of these Lagragians are a se- 
ries of terms with increasing chiral dimension, 
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where the ellipsis stand for terms of yet higher order. The 
explicit expressions for the various terms are well docu- 
mented in the literature, see e.g. [14]. The transition op- 
erators derived from this effective Lagrangian are then 
sandwiched between wave functions that were consistently 
generated from chiral nuclear effective field theory. The 
latter we take from the recent work of Ref. [15]. The vari- 
ous contributions can be organized in terms of a consistent 
power counting in terms of a small parameter q, like e.g. 
a meson mass, energy or nucleon three-momentum (with 
respect to the typical hadronic scale, say the mass of the 
rho meson). For a generic matrix element one has, 

M = q v F{llti, (5) 

where /i is a renormalization scale and the function T 
is of order one. Furthermore, v is a counting index, i.e. 
the chiral dimension of any Feynman graph, for the case 
of pion production off nuclei, see e.g. [16]. In terms of 
this counting index, we include all terms with v = and 
v = — 1 (see also [3] for details). 

The transition matrix elements are generated by two 
very different types of contributions, 

M = M ss + M th , (6) 

where "ss" and "tb" denote the single-scattering and the 
three-body contribution, respectively. Here, single scatter- 
ing means that the pion is emitted from the same nucleon 
to which the photon couples with the other nucleon act- 
ing as a mere spectator. Processes involving both nucle- 
ons are often called exchange currents, but we follow here 
the notation due to Weinberg [17]. At third order, all tb 
diagrams involve graphs with one pion in flight, whereas 
at the order considered here, four-nucleon contact terms 
with the photon absorption and pion emission from dif- 
ferent nucleon legs come in (as discussed in more detail 
below). These terms can be understood from integrating 
out heavier mesons, see e.g. [18,19]. In such a picture, such 
diagrams would thus correspond to heavy meson exchange 
currents. 

The ss terms are of course sensitive to the elementary 
proton and neutron electroproduction multipoles, prop- 
erly boosted to the pion-deuteron center-of-mass frame 



(for details, see [3]). The proton amplitude is fixed from 
the chiral perturbation theory study of Ref. [10] . As in our 
earlier work, we perform two types of fits, where we have 
one respectively two undetermined parameters related to 
the elementary nir° amplitude. In the fits of type 1, we 
have one free fifth order parameter, called a", and we 
include the constraint from a leading order low-energy 
theorem to the fourth order counterterms (as explained 
in [10]). The numerical value for the LEC from reso- 
nance saturation is -4.11 GeV~ 4 (comprised of the con- 
tribution from the Z\(1232) in static approximation and 
vector mesons). For the fits of type 2, we relax this con- 
straint and thus have two fit parameters, the low-energy 
constants and a". These two fitting procedures give us 
a measure of the theoretical uncertainty at the order we 
are considering. We also get another measure of the the- 
oretical uncertainty by considering chiral EFT at NNLO 
and varying the cut-off in the Lippmann-Schwinger equa- 
tion. This will lead to bands for the various observables 
rather than to lines as e.g. in [3]. We will come back to 
this topic when discussing the results. We note already 
here that the uncertainty related to the two fit procedures 
is sizeably bigger than the one induced by the cut-off vari- 
ation. 

In what follows, we will use as kinematical quantities 
the virtuality of the photon, k 2 , which is negative in elec- 
tron scattering (in the literature one often uses the pos- 
itive quantity Q 2 — —k 2 ), the photon polarization e and 
the pion excess energy, AW, that is the energy of the pro- 
duced pion above threshold in the ir — d center-of-mass 
system. For a more detailed discussion of the kinematics, 
see e.g. [1]. 

2.2 Fourth-order three-body contributions 

At third order, there are 8 tree graphs contributing to 
the three-body corrections, see Fig. 4 of Ref. [3] (using 
Coulomb gauge). At fourth order, there are altogether 59 
non- vanishing diagrams, as shown in Fig. ^ These are 
tree graphs with exactly one insertion from the dimension 
two chiral pion-nucleon Lagrangian, C n L. More precisely, 
in that figure we have only shown the topologically in- 
equivalent diagrams, the numbers under certain graphs 
denote how many diagrams can be generated if one at- 
taches the pion emission vertex on the left or the right 
nucleon line above or below the pion exchange line. There 
is also a whole new class of short-distance diagrams includ- 
ing the leading (momentum-independent) four-nucleon in- 
teractions (the last four diagrams in Fig. ^ an d the cor- 
responding Okubo-corrections, i.e. the diagrams with two 
close- by energy denominators, are not shown. For details 
on this point, see [20]. A more general discussion one 
these re-orthonormalization diagrams can be found e.g. 
in [21,22].). The corresponding four-nucleon LECs have 
already been determined in the fits to nucleon-nucleon 
scattering data. We use here the values collected in ta- 
ble ^ for the corresponding NNLO wave functions. 
However, these values stem from a fit to the NN phase 
shifts at low energies. A fit including pion production data 
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LEC 


A = 450 MeV 


A = 650 MeV 


Ciso [lO" 4 GeV" 2 ] 
Casi [10" 4 GeV- 2 ] 


-0.151 
-0.168 


-0.149 
-0.130 



Table 1. Values of the four-nucleon LECs for the NNLO wave 
functions of [15] used here. A is the cut-off in the Lippman- 
Schwinger equation (the spectral function cut-off is fixed at 
A = 650 MeV) 



(which is from the kinematical point of view closer to the 
process considered here) can lead to an increased theo- 
retical uncertainty in the determination of these LECs. 
To make that point more transparent, we briefly discuss 
the determination of the leading ANir LEC related to the 
D-term in [25]. From low-energy observables one obtains 
for this coupling in dimensionless units, Cd/(F%A^) = 
1.8 . . .3.6, whereas its determination from P-waves in the 
reaction pp — > ppir gives somewhat lower (but still consis- 
tent) values, cd = 1.4. . .2.1 [26] 1 . Also, at fourth order 
we will have additional P-wave LECs from the single nu- 
cleon amplitudes. With these two observations in mind, 
we will also perform fits were we leave Ciso and C^si 
as free parameters. This will be discussed in more de- 

(2) 

tail when we present the results. Furthermore, C^^ con- 
tains some terms with fixed couplings (due to the con- 
straints of Lorentz invariance) and other terms with fi- 
nite and scale-independent low-energy constants, denoted 
Cj. These LECs can e.g. be determined from the analysis 
of low-energy elastic pion-nucleon scattering data utiliz- 
ing chiral perturbation theory. From the hadronic cou- 
plings only the term ~ C4 contributes, we use the value 
c 4 = 3.4 GeV- 1 [23]. 

We note that all the diagrams shown in Fig. ^ are 
irreducible, where we use the following definition for ir- 
reducible diagrams: A diagram is called irreducible, if it 
contains no contributions from the two-nucleon potential 
(see e.g. the related discussion in [24]). This definitions 
becomes obvious if one considers the Lippmann-Schwinger 
equation, 

(6,7r°|T|& 7 *) = 
(<j)d 7T° |(K-jv + V„nn) + (Vkn + Vknn) 

1 

X — — — — — — — — 

E — H — Vnn — Vkn ~ 'ttNN — y~iN — VyNN + *e 

x(V-, N + V 7NN )\c/) d -f*} ,(7) 

and the deuteron wave function <fi d is obtained from the 
solution of the Schrodinger equation, 



(H + V NN ) \4> d ) = E d \(f> d ) 



(8) 



Here, V^n, VyN, ■ ■ ■ are the pertinent pion-nucleon, photon- 
nucleon, . . . transition potentials subject to the chiral ex- 
pansion as explained above. This concept of reducibility is 



1 Note that we have a different convention for the sign of the 
axial- vector coupling The relation between the coupling 
constant 8 used in that paper and cd is given by 8 — — 0.143cr>. 



depicted in Fig. [3 where two typical reducible diagrams 
are shown, which include either the lowest order one-pion- 
exchange or the short-distance parts of the two-nucleon 
interaction. All diagrams depicted in Fig. ^ fall essen- 
tially in four classes, labeled a), b),c) and d). The first 
two classes were already present at third order [3], these 
are the seagull- type and the pion- in-flight type graphs, 
respectively. Class c) collects the so-called time-ordered 
graphs, were the photon couples to a nucleon line and the 
pion is emitted from a nucleon while the exchanged pion is 
in flight. All diagrams build from the four-nucleon contact 
terms are collected in class d) . Representative diagrams for 
these four classes are shown in Fig.|3 In fact, a) falls into 
two subclasses, depending on the value of the energy q' of 
the exchanged pion. One either has q' = + 0(1/ m) or 
l'o = 1o+0(l/m), with go the energy of the produced pion. 
The diagrams are evaluated using Fourier-transformation 
techniques, as briefly discussed in Add.IAI We also remark 
that all our coordinate space integrals are finite due to the 
exponential fall-off of the chiral EFT wave functions. In 
contrast to what was done in [4], we thus do not need 
to introduce an additional cut-off in the fourth order tb 
terms. 



3 Results and discussion 

In this section, we display the results for the multipoles, 
differential and total cross sections and the S-wave cross 
section ao d for the two fit strategies. We have performed 
calculations with the chiral EFT wave functions at NNLO 
[15] for cut-offs in the range from 450 to 650 MeV (with the 
spectral function cut-off A fixed at 650 MeV, for details 
see [15]). Consequently, for the observables we will obtain 
bands rather than single line. This method of estimating 
the theoretical uncertainty is one of the advantages of the 
chiral effective field theory approach employed here. We 
note that at threshold we have performed the calculations 
in momentum and in coordinate space, which serves as a 
good check on the nontrivial numerical evaluation of the 
various integrals. 



3.1 Results with fixed four-nucleon LECs 

We first discuss the result obtained keeping the four-nucleon 
LECs Ciso and C3S1 fixed at the values given in sec- 
tion In table [21 we collect the fitted LECs related to 
the elementary neutron amplitude. As discussed before, 
we have performed two types of fits. We see from the ta- 
ble that the wave function dependence of the LECs is very 
weak, which points towards the conclusion that the perti- 
nent matrix elements are dominated by the contributions 
from the long-ranged pion exchange. The resulting values 
are also not very different from the ones obtained in [3], 
which are also given in the table. 

Before showing results for the pertinent observables, we 
discuss a few general features of the new contributions. 
First, we observe that the 4N contact interactions con- 
tribute mostly to the P-waves, their S-wave contribution 
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LECs A = 450 MeV A = 650 MeV 0(<j 3 ) 

ag [GeV -4 ] 4.148 4.140 4.767 

al [GeV -4 ] -7.085 -6.730 -5.644 

aV; [GeV~ 5 ] -38.71 -37.05 -27.51 

Table 2. Values of the fitted LECs for the NNLO wave func- 
tions from [15]. The values for refer to the fits 2, whereas 
the corresponding as belongs to the respective fits 1. For com- 
parison, the LECs from [3] are also shown. 



is very small. Further, these P-wave contributions mostly 
feed into the magnetic M-Jq multipole and cancel to some 
extent the corresponding contribution from the third or- 
der tb terms. Second, the other new diagrams also lead 
to more changes in the P-wave multipoles, simply since 
at this order there are no free parameters. We also note 
that the value of the constant in the fits 2 agrees nicely 
in magnitude with the resonance saturation estimate, but 
comes out with opposite sign. 

In Figs. 14151 we show the differential cross sections for 
fits 1 and 2 employing the NNLO wave functions in com- 
parison to the MAMI data [1]. These two bands (which 
are generated by utilizing the NNLO wave functions with 
the cut-off A = 450 MeV and 650 MeV) corresponding to 
the two fit procedures can be considered as a measure of 
the theoretical uncertainty at this order. We note that the 
bands due to the wave function dependence are very thin, 
which shows that the effective field theory calculation of 
these wave functions is of sufficient accuracy. The uncer- 
tainty generated from the two fit procedures is comparable 
to the experimental errors. We remark that the differences 
between the two fit procedures is somewhat smaller as 
it was the case when only the third order tb corrections 
were included [3], see the dashed lines in Figs. 14151 Con- 
sequently, as it should be in a converging effective field 
theory, the theoretical uncertainty has become smaller as 
compared to Ref. [3], although this improvement is mod- 
erate. 

The corresponding total cross sections as a function 
of the excess energy AW and of the photon polarization 
e are shown in Figs. |5] for fit 1 and in Figs. for fit 2 
and the NNLO wave function with A ranging from 450 
to 650 MeV, respectively. We notice that for fit 1 with 
increasing excess energy and, in particular, with increasing 
photon polarization the data are systematically below the 
chiral prediction. Due to the fitting procedure, the slopes 
of the various curves for the Roscnbluth separation shown 
in the right panel of Fig. are of course correct. These 
results are very similar to the ones obtained in [3] for 
fit 1. We note that the predictions for fit procedure 2 are 
visibly improved compared to the third order calculation. 
This is further seen by looking at the transverse and the 
longitudinal threshold S-wave multipole shown in Fig. |H1 
The prediction for \Ld\ for fit 2 is now within the error 
bar of the MAMI result (which was not the case in [3] ) . 
We note that the real part of the predicted longitudinal 
amplitude is negative, consistent with our findings in [3]. 
The resulting S-wave cross section aod is shown in Fig. 



The theoretical uncertainty is a bit smaller than in case 
of the third order tb calculation [3] . 

3.2 Results with fitted four-nucleon LECs 

As argued before, we will also perform fits leaving the 
values of the two four-nucleon LECs as free parameters. 
Since as we noted before these operators contribute very 
little to the S-waves, and not at all at threshold, the fits of 
type 1 are performed in two stages. First, a§ is adjusted 
to give the proper value of |L<j| at k 2 = —0.1 GeV 2 . Then, 
the two four-nucleon LECs are fitted to the total cross 
sections. The fits of type 2 are performed as before, all pa- 
rameters are fixed on the total cross sections. Having said 
this, we collect in table El the values of the various LECs 
for the two types of fits. We note that the four-nucleon 
LECs come out much larger than from the investigation 
of nucleon-nucleon scattering, cf. tabled However, as we 
discussed earlier, leaving these LECs as free parameters 
effectively subsumes some effects from the single nucleon 
P-wave contributions, so that these fits should be consid- 
ered indicative only. 



LECs 


A = 450 MeV 


A = 650 MeV 


a? [GeV" 4 ] 


4.500 


4.751 


al [GeV" 4 ] 


-8.549 


-8.573 


Ciso [lO" 4 GeV" 2 ] 


-0.589 


-0.664 


C 3S i [lO" 4 GeV -2 ] 


-1.168 


-1.105 


a? [GeV- 5 ] 


-38.71 


-37.05 


Ciso [10- 4 GeV" 2 ] 


-0.325 


-0.387 


C 3 si [10" 4 GeV- 2 ] 


-0.762 


-0.890 



Table 3. Values of the fitted LECs for the NNLO wave func- 
tions from [15]. The values for 03,4 refer to the fits 2 with the 
corresponding values for Ciso and C3S1, whereas the corre- 
sponding ag belongs to the respective fits 1, again with the 
corresponding four-nucleon LECs. 



The resulting total cross sections are shown in Fig. ITU1 
We observe a clear improvement for fit 1 and a very good 
description for fit 2. We remark, however, that the corre- 
sponding differential cross sections for AW > 1.5 MeV are 
more symmetric around cos((9) = (bell-shaped) than the 
ones shown in the preceding section, whereas the data in- 
dicate a peaking into the backward direction. This points 
towards an insufficient accuracy in the description of some 
of the P-wave multipoles. The resulting longitudinal deuteron 
S-wave multipole is shown in Fig. El While \Ed\ is almost 
unaffected, \L d (k 2 = -0.1GeV 2 )| for fit 2 agrees with the 
empirical value. Consequently, the band for the S-wave 
cross section is much narrower than before, cf. the right 
panel of Fig. ^2 This is because the cross sections are 
much sensitive to the longitudinal S-wave amplitude due 
to the kinematical enhancement proportional to the lon- 
gitudinal polarization el ~ 9, see [1]. 
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4 Summary and outlook 

In this manuscript, we have considered neutral pion elec- 
troproduction off deuterium in the framework of heavy 
baryon chiral perturbation theory, extending and improv- 
ing upon the work presented in [3]. The salient results of 
this study are: 

i) We have calculated the fourth order three-body cor- 
rections. These consist of 51 one-pion exchange dia- 
grams with exactly one insertion from the dimension 
two pion-nucleon Lagrangian and 8 diagrams with one 
insertion from the lowest order four-nucleon interac- 
tion Lagrangian. In principle, all parameters are fixed 
from earlier studies of pion-nucleon and nucleon-nucleon 
scattering. The deuteron wave functions are taken con- 
sistently from the recent chiral EFT study of [15]. 

ii) As in the earlier work [3] , in which the three-body cor- 
rections were only considered at third order, we have 
performed two types of fits to the MAMI data [1] . The 
results are similar to the ones found there, although 
we observe a moderate improvement in the theoreti- 
cal uncertainty. This is most pronounced for the lon- 
gitudinal S-wave multipole \Ld\ in fit 2. We have also 
considered the dependence on the deuteron wave func- 
tions related to the cut-off in the Lippmann-Schwingcr 
equation. It turned out that this dependence is very 
weak. Therefore, the process is dominated by long- 
range pion physics and thus sensitive to the elementary 
7r°n amplitude. 

iii) We have also performed fits where we have left the 
four-nucleon LECs as free parameters, based on the 
argument that in that way one can effectively subsume 
new P-wave LECs from the single scattering contribu- 
tion. This leads to a visibly improved description of 
the total cross sections and the longitudinal S-wave 
multipole. However, the differential cross sections come 
out bell-shaped, in contrast to the experimental find- 
ings [1]. 

In conclusion, we have further sharpened the theoreti- 
cal framework to analyse neutral pion production off deu- 
terium in a model independent framework. It remains to 
be seen whether a further improved description of the data 
will be possible when the fourth order pion-nucleon P- 
wave multipoles at fourth order become available. 
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A Fourier transformations 

The diagrams are evaluated in momentum space. To cal- 
culate the pertinent matrix elements, we perform Fourier 



integrations, see e.g. [4]. The typical structures to be eval- 
uated take the form 

p' vl ...p' vm ^ ( P >) 

= (2^72 / dVe-*'-* 1 (-*ft,i)...(-i&™)0(ri) 

Pal . . .p ak $(p) 

where p (p') are initial (final) state relative nucleon mo- 
menta, <p is the coordinate space deuteron wave function 
and similarly <P denotes the momentum space wave func- 
tion. As stated, we have to consider 5 classes of diagrams, 

1) Diagrams of type a) with q' = q n + 0(l/m), 

2) Diagrams of type a) with q' = 0(l/m), 

3) Diagrams of type b), 

4) Time-ordered diagrams of type c), 

5) The NN-diagrams of type d). 

In the following, symbolizes an arbitrary spin-operator. 

Consider first the diagrams 1). The momentum of the 
exchanged pion is 



. . k q 

q >=p-pi + - + -. 



(A.2) 



From the momentum space expression for this class of 
diagrams, one deduces the following pertinent Fourier- 
transform 



1 



q' ■ q' + u)% - <7q - it 4n 
with 



in J 



-S r 



drre 



(A.3) 



(A.4) 



Here, lo c is the physical value of the nmr + and the ppir~ 
threshold, respectively. So then we have to evaluate the 
following integrals 



J d 3 pd 3 p'p' ul ...p' um <P* (P')^7 



1 



q' ■ q' + w% - ql - it 







xp al ...p ak $(p) = 2tt 2 J d 3 r(-id ul )...(-id um )(f)(r) 



-S r 



id <7l ...id ak cf>(r) e~ 4 (*+*)- 



(A.5) 



Next, we evaluate the diagrams 2). The momentum of 
the exchanged pion now is 



k 

2 



q =p p 

and the pertinent Fourier integral is 
1 



(A.6) 



■ q' + Ml + - it 



1 r , o~ M „+ r 
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leading to the structures 
xp a x...p ak ${p) = 2?r 2 / Sri-idyi) . . .{-id vm )4>{r) 



X0- 



-id a i . . . id ak (j> (r) e~ 



«(*-*)• 



(A. 



We now consider the diagrams 3). Here, the photon 
couples to the pion in flight, thus we have two propagators. 
The momentum of the exchanged pion before the photon 
absorption is 



, . k q 

q'=P~P +2 + 2 



(A.9) 



and after the photon absorption this momentum is changed 
to 

q" = q'-k. (A.10) 
We have two types of Fourier-transforms, these are: 

1 

(q" ■ q" + lo 2 - ie)[q' -q' + to^-q 2 - ie) 

rl l 

dz - , (A.ll) 



" (I 2 + m' 2 - ie) 2 



with 



and 



I = q — zk , 

+ (l-z) (5 2 +k 2 z) , (A.12) 



/2 2 
Til = 0t> Z 



Qfll l'fj.2 



[q" ■ q" + lu 2 - ie){q' -q' + u 2 -q 2 - ie) 



= / dz 

/g Sum' 



d 3 re- llr 



x [—id pa + zk^)[-id^ + zk^) < 



(A.13) 



Note that this function is still square integrable. Higher 
order derivatives, however, can not be treated in this fash- 
ion. The corresponding integrals are thus 



J d 3 pd 3 p'p' lyl ...p' vm <t>* [p') 



[q" ■ q" + uj 2 - ie)[q' ■ q' + iv 2 - q 2 - ie) 

f 1 7T 2 f 
x0p al . . .p ak ${p) = / dz— d 3 r 

Jo m J 
x [-idvi) ■ ■ ■ {-idum)^ [r) [-id^i + zk^) 

x(-id^ 2 + zk^) c~ m r idcri . . . id ak 4> ( r ) 

xe -Ki~( z ~T;) k )- r . (A.14) 

We now discuss the diagrams 4), which have the momen- 
tum transfer given either by Eq. I|A.2(I or by Eq. I|A.6(I . 



Due to the cancellation with the corresponding Okubo- 
type diagrams, no new structures as compared to the other 
type of diagrams appears (for a more detailed discussion 
on this point see [20]) 

Finally, we turn to the NN diagrams 5). Since these only 
involve contact interactions, we have to consider simple 
structures of the form 

J d 3 p d 3 p' p' vl . . . p' vm $* (p') 0p aX . . . Pak <P (p ) 

= (2tt) 3 (-idui) . . . [-id vm )4> (n) | n=0 

x (^...(^^(rj)^ (A.15) 

We end this appendix with a note on how we treat 
the contributions arising from the exchange of nucleon 1 
with nucleon 2 when calculating the pertinent amplitude, 
A = A[l, 2) + (1 <-» 2). The corresponding momenta of the 
nuclcons in terms of the center-of-mass and the photon 
and pion momenta are 



Px=P 



k 

2 



Pi=P 



2 



P2 = ~P 



k 

2 



p' 2 = -P' 



2 



(A.16) 



Under a parity transformation p — > p and p 1 — » — p', 
the nucleon momenta behave as p\ — > p 2 , J>2 —> Pi , 
P'l P 2 , P2 Pi > which gives 

J d 3 pd 3 p'<P* (p') / [at , <r 2 , p x , p 2 , p[, p' 2 ) $ (p) 

= J d 3 pd 3 P '^ ( P ')f( ( T 1 , ( T 2 ,p 2 ,p 1 ,p' 2 ,p' 1 )<P(p). 

(A.17) 

Thus, matrix elements are invariant under operations of 
the form 

f ((T 1 ,(T 2 ,p 1 ,p 2 ,p' 1 ,p' 2 ) -» 

/(<Ti,CT 2 ,P2,Pl,P2,Pi) • 

Consequently, the operation (1 <-> 2) need only be applied 
to the spin and isospin operators and leave all other quan- 
tities unchanged. Furthermore, because of isospin symme- 
try we only have the operators t% ■ t 2 — t^t 2 and rfrf, 
which are symmetric under the interchange (1 <-> 2). There- 
fore, we finally only need to change the spin indices when 
calculating the contributions from the exchange diagrams. 



B Angular integrations 

Here, we briefly discuss how to efficiently perform the 
angular integrations of the Fourier integrals discussed in 
App-El We define the following matrices constructed from 
the momentum space deuteron wave functions, 



in,m[r) — 



1 



dqq 2 u[q) [qr) n j m [qr) , 
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n 1/2 

w n m(r) = : — 

''»"*\ ' r n+m 



dqq 2 w{q) {qr) n j m {qr) , 



(B.l) 



where u(q) and w(q) are the momentum-space S- and D- 
wave components of the deuteron wave function, and the 
in (qr) are spherical Bessel functions. The usual coordinate 
space expressions can be obtained using matrix elements 
u ,o and Wq,2, i-e. 



uq,o 



u(r) 



W0.2 



w(r) 



(B.2) 



This notations is particularly useful since we have to deal 
with derivatives of the wave functions. Using the prop- 
erties of the spherical Bessel functions, one obtains the 
following recursion relation 

r 2 u n+ i. m+2 (r) - (2m + 3)M„ >m+ i(r) + u n+ i, m (r) = . 

(B.3) 

This allows one to express the derivatives of uq^q in the 
following way 

u o,o0) = -ru hl (r) , 

u o,o( r ) = -«l,l( r ) +i' 2u 2.2( r ) , 

<; (r) = 3ru 2 , 2 (r) - r 3 u 3 , 3 (r) . (B.4) 

Similarly, the derivatives of iuo,2 can be expressed as 

w 'o,2( r ) = ~- (5w Qt2 (r) -toi,i(r)) , 

w o,2( r ) = "J (3Qw ,<i(r) - 6wi,i(r) - r 2 io 2 2(r)) 

< 2 (r) = i (-2Khuo,a(r) +42iBi,i(r) 

+6r 2 w 2 , 2 (r)+r 4 u;3,3(r-)) . (B.5) 

Note that one only has to calculate diagonal matrix ele- 
ments. 
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Fig. 1. Fourth order diagrams for the three-body corrections (in Coulomb gauge). Solid, 
dashed and wiggly lines denote nucleons, pions and photon, in order. The heavy dot stands 
for an insertion from the dimension two pion-nucleon Lagrangian. The number under certain 
graphs stands for the number of topological equivalent graphs as explained in the text. The 
ellipsis at the end stands for Okubo-type corrections to the graphs involving the four-nucleon 
contact interactions and also one-pion exchange. 
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Fig. 2. Classes of reducible diagrams as explained in the 
text. The shaded area denotes the two-nucleon interaction, 
which includes either the leading one-pion exchange (left) of the 
momentum-independent four-nucleon interactions (right). Such 
types of diagrams are consistently not included in the three-body 
interaction kernel. 




a) b) c) d) 



Fig. 3. Classes of irreducible diagrams as explained in the text, a) represents a seagull, b) 
a pion-in-flight, c) a time-ordered and d) a short-distance diagram, in order. 
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Fig. 4. Differential cross section at AW = 0.5 MeV (left column) and AW = 1.5 MeV (right column) at three 
different values of the photon polarization for the NNLO wave functions in comparison to the MAMI data [1]. 
Fits 1 and 2 are given by the blue (upper) and the red (lower) bands, respectively. The bands are generated by 
varying A from 450 to 650 MeV. The dashed lines are the corresponding results from [3]. 
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Fig. 5. Differential cross section at AW = 2.5 MeV (left column) and AW = 3.5 MeV (right column) at three 
different values of the photon polarization for the NNLO wave functions in comparison to the MAMI data [1]. 
Fits 1 and 2 are given by the blue (upper) and the red (lower) bands, respectively. The dashed lines are the 
corresponding results from [3]. 



H. Krebs et al.: Improved analysis of neutral pion electroproduction off deuterium in chiral perturbation theory 




1 2 3 

AW [MeV] 



to 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

8 



Fig. 6. Left panel: Total cross section as a function of AW for three different 
values of the photon polarization in comparison to the MAMI data [1] for fit 1 
and the NNLO wave functions. The upper/middle/lower band corresponds to the 
largest /medium/smallest value of e. Right panel: Total cross section as a function 
of the photon polarization e for four different values of the pion excess energy AW . 
The highest band corresponds to the largest value of AW and correspondingly for 
the other bands/values of AW . 




Fig. 7. Left panel: Total cross section as a function of AW for three different 
values of the photon polarization in comparison to the MAMI data [1] for fit 2 
and the NNLO wave functions. The upper/middle/lower band corresponds to the 
largest /medium/smallest value of e. Right panel: Total cross section as a function 
of the photon polarization e for four different values of the pion excess energy AW . 
The highest band corresponds to the largest value of AW and correspondingly for 
the other bands/values of AW . 
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Fig. 8. Threshold multipoles \E&\ (left panel) and \Ld\ (right panel) as a function of 
the photon virtuality in comparison to the photon point data from SAL [27] and the 
electroproduction data from MAMI [1]. The sign of the experimental result for Ld is 
taken to agree with the theoretical prediction. Ed'- The blue (upper)/red (lower) band 
refers to fit 1,2, in order (for Ld, the upper/lower refers to the value at Q 2 = 0.1 GeV 2 ). 
The dashed lines are the corresponding results from [3]. 
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Fig. 9. The S-wave cross section a^d as a function of the photon virtu- 
ality. Solid bands (dashed lines): Range obtained with the fourth (third) 
order three-body corrections. The dashed lines are the corresponding 
results from [3]. 
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Fig. 10. Left panel: Total cross section as a function of AW for three different 
values of the photon polarization in comparison to the MAMI data [1] for fit 1 
and the NNLO wave functions. The upper/middle/lower band corresponds to the 
largest /medium/smallest value of e. Right panel: Same for fit 2. 




Fig. 11. Fits with varying four-nucleon LECs. Left panel: Threshold multipole \Ld\ 
as a function of the photon virtuality in comparison to the electroproduction data 
from MAMI [1] . The sign of the experimental result for Ld is taken to agree with the 
theoretical prediction. The blue (upper)/red (lower) band refers to fit 1,2, in order. 
Right panel: The S-wave cross section aod as a function of the photon virtuality. The 
solid bands give the range obtained with the fourth order three-body corrections, 
dashed lines show the third order tb results of [3]. 



